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ABSTRACT 

We present the search methods and initial results for transient radio signals in the High 
Time Resolution Universe (HTRU) Survey. The HTRU survey's single-pulse search, 
the software designed to perform the search, and a determination of the HTRU survey's 
sensitivity to single pulses are described. Initial processing of a small fraction of the 
survey has produced 11 discoveries, all of which are sparsely-emitting neutron stars, 
as well as provided confirmation of two previously unconfirmed neutron stars. Most 
of the newly discovered objects lie in regions surveyed previously, indicating both the 
improved sensitivity of the HTRU survey observing system and the dynamic nature 
of the radio sky. The cycles of active and null states in nulling pulsars, rotating radio 
transients (RRATs), and long-term intermittent pulsars are explored in the context 
of determining the relationship between these populations, and of the sensitivity of 
a search to the various radio-intermittent neutron star populations. This analysis 
supports the case that many RRATs are in fact high-null-fraction pulsars (i. e. with null 
fraction >0.95), and indicates that intermittent pulsars appear distinct from nulling 
pulsars in their activity cycle timescales. We find that in the measured population, 
there is a deficit of pulsars with typical emission timescales greater than ^300 s that 
is not readily explained by selection effects. The HTRU low-latitude survey will be 
capable of addressing whether this deficit is physical. We predict that the HTRU 
survey will explore pulsars with a broad range of nulling fractions (up to and beyond 
0.999), and at its completion is likely to increase the currently known RRATs by a 
factor of more than two. 
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1 INTRODUCTION 

Radio single-pulse searches within the last decade have re- 
vealed a number of fascinating objects within and outside 
of our own Galaxy. Pulsed radio emission originating from 
other galaxies could give an estimate of the combined free 
electron content of their host galaxies and of the intergalac- 
tic medium; thus far, searches have detected single pulses 
from neutron stars in two close satellites of our Galaxy, 
including gi ant pulses from B0540-69 in the Large Magel- 
lanic Cloud ((Johnston fc Romanil2003l ) and a probable tran- 
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sient neutron star in the Sculptor Spheroidal Dwarf Galaxy 
(|Rubio-Herreral [2OI0I ) . A number of other nearby galaxies 
have been the target of single-pulse searches, showing only 
sparsely detected events of inc onclusive origin (e. g. ga l axies 
M33 and M31, searched by iMcLaughlin fc Corded |2003| 
and Rubio-Herrera 2010l . respectively). The discovery by 
iLorimer et al.l ( 2007 ) of a seemingly dispersed, 30- Jy im- 
pulse (at a dispersion measure, DM, of 375pccm~^) ap- 
peared to represent the first discovery of an extragalac- 
tic pulse with high significance that was of non-neutron 
star origin. However, the recent discovery of pulses of am- 
biguous terrestrial origin with frequency-sweeps that mimic 
the cold plasma dispersion relation and primarily appear 
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around DM ~ 370pccm~'^, have cast some doubt on the 
Lorimer et al. pulse as being extragalactic and provided 
an a dditional, terrestrial target for single-pulse searches 
jBurke-Spolaor et al. 2011, '). 

The most frequently detected sources of transient radio 
emission are pulsars. Single-pulse searches in recent years 
uncovered what was la belled the "rotating rad io transient" 
(RRAT) phenomenon ()McLaughhn et al.1 '20061. which en- 
capsulates sparsely radio-emitting neutron stars, that due to 
their sporadic emission, were undetectable by the Fourier- 
based search techniques currently in standard use for pul- 
sar surveys. The point has recently been made that in fact 
the term RRAT is rather a "detection label" than a de- 
scription of a physical phenomenon: Keane (2010a) notes 
that whether a discovery is considered "RRAT" is highly 
dependent on both survey length and an object's rota- 
tional period. This and other mounting evidence indicate 
that RRATs are unlikely to have a physically distinct ori- 
gin from radio pulsars. This includes, for instance, that 
RRAT pulse energy distributions tend to obey log-normal 
distributi ons similar to bot h non-nulling and nulling pul- 
sars (e. g. ICairns et aLlliooH . Keane 2010b; Miller et al. sub- 
mitted) and that there is a lack of distinction between 
RRAT and p ulsar Galactic distributions an d pulse width 
distributions l|Burke-Spolaor fc Bailed l2010t ) . Additionally, 
if RRATs were a distinct phenomenon unrelated to other 
neutron star populations in an evolutionary sequence, the 
implied birthrate of neutron stars would far exceed the rat e 
of supernovae that produce them (|Keane fc Krameill2008f ). 
As such, the current understanding of the objects discov- 
ered preferentially in single-pulse searches is that they are 
a mix of both mod ulated pulsars with lon g-tailed pulse 
energy distributions (jWeltevrede et al] |2006| ) , and pulsars 
at t he m ost extreme end o f the nulling pulsar popula- 
tion iBurkc-Spolaor & Bailel l2010l : Keane 2010b; Miller et 
al. submitted). 

However, questions remain: why do the objects appear 
to have a period-derivative and/or magnetic field distribu- 
tion that si ts higher than a verage pulsars of the same pe- 
riod range (|McLaughlin et a l. 2009)? What are the statis- 
tics of peculiar phenomena in these extreme nulling ob- 
jects such as the RRAT-pulsar mode switching of PSR 
J0941-39 l|Burke-Spolaor fc Bailed l2010t ). the multi-modal 
(and latitude-dependen t) behaviour of PSR J1119-6127 



Weltevrede et al.1 I2OI0I), or the glitch activity of PSR 



J1819-1458 l|Lvne et al.l |2009| )? What causes nulling and 



what is the distribution of pulsar nulling fractions (again 
holding implications for the neutron star birthrate)? And 
finally, do (and if so, how do) the most extreme nulling 
pulsars fit into an evolutionary progression between av- 
erage radio pulsars, nulling pulsars, ra dio quiet neutron 
stars , and the magnctar po pulation (e. g. [M cLaughl in et au 
l2009l : lLvne"et al. 2009; Bur ke-Spolaor fc B ailcs 201(i; Keane 
2010b)? Acquiring larger statistical sample of these ob- 
jects is among the next essential steps in understand- 
ing RRATs and transient radio neutron star phenom- 
ena (currently a tot al of ~40 such objects are known, 
Hcssels et al. 2008: Deneva et al. 2009; Keane et al. 2010; 
Burke-Spolaor fc BaiQl2010l : lRubio-Herrer3l201oi ). 



tions 5 <+lQ° l|Keith et al.1 l20ld : hereafter HTRU Pa- 
per 1). The survey employs a new digital backend, the 
"Berkeley- Parkes-Swinburne Recorder" (hereafter BPSR), 
that has been installed at the Parkes Telescope for the 20- 
cm multibeam receiver l|Stavelev-Smith et al.lTl99^ . BPSR 
allows improved digitisation levels, frequency and time res- 
olution over the previous analogue instrument that has 
been used for previous southern single-pulse search es and 
studies (i. e.lBurke-Spolaqr fc Bailes 2010 . McLaughlin et al.l 
(|2006l)lKeane et al.1 (|2010l )). This affords the HTRU survey 
unprecedented sensitivity to sub-second duration, dispersed 
single impulses of radio emission in the southern sky. This 
paper describes the techniques used in the survey to search 
for single pulses. 

As an all-sky survey, the HTRU survey is sensitive 
to terrestrial. Galactic, and extragalactic sources of radio 
pulses. The HTRU intermediate and high latitude surveys 
will cover Galactic longitudes —120° < I < 30° and latitudes 
|6| < 15°, and declinations S < +10° (not covered by inter- 
mediate pointings) , respectively. Although no processing has 
yet commenced for the HTRU low-latitude survey, that sur- 
vey will cover areas \h\ < 3.5° and —80° < I < 30°, with 
70 minutes per pointing. A parallel northern effort is under- 
way at Effelsberg Radio Telescope, affording a full sky sur- 
vey for pulsars and transients. The low and intermediate lat- 
itude surveys share common areas of sky coverage to previ- 
ous surveys which employed the analogue filterbank installed 
on the Parkes 20-c m Multibeam receiver. The "Parkes 
Muhibeam Survey" l|Manchester et all [2OO1I ') covered |fe| < 
5° with 35-minute pointings, and was searched for singl e 
pulses bv lMcLaughlin et all (|2006h and lKeane et al.l |20l3), 
revealing the bulk of the currently recognised RRATs. Here- 
after we refer to those two searches collectively as the 
PKSMB searches. At higher Galac tic latitudes, a single- 
pul se search was perfo rmed on th e lEdwards et ahi (l200ll) 
and Jacobv et al.l (|2009l ') surveys bv lBurke-Spolaor fc Bailea 



(hereafter BSB). Despite the common sky coverage 
and high success rate of the BSB and PKSMB searches, 
we nevertheless expect to discover new sources of pulsed 
radio emission in these areas; some pulsing radio sources 
may not have been previously detected due to the source's 
low event rate, insufficient sensitivity in previous surveys, or 
the improved dynamic range of the BPSR observing system 
over the previous analogue filterbank system installed on the 
Parkes multibeam receiver, whose dynamic range was hurt 
by \/]V-statistics (explored in more detail in i]2. 21 below). 



The High Time Resolution Universe (HTRU) Sur- 
vey is the first digital, all-southern-sky survey for pul- 
sars and fast (sub-second) transients, covering declina- 



In this manuscript, we present the search methods, 
sensitivity, and current status of the HTRU survey single- 
pulse search (321 • Because this survey required concurrent 
searches using periodicity and single-pulse techniques, it af- 
fords the opportunity to explore efficient modes of search 
operation, attempting to perform these techniques in par- 
allel with a minimum of duplicated computing effort. We 
also report on initial discoveries in the intermediate- and 
high-latitude surveys (331) • In SI we explore the detectabil- 
ity of our discoveries in previous surveys of overlapping re- 
gions, present our new discoveries in terms of the window- 
ing description of pulsar intermittence presented in BSB, 
and make predictions for the sparsely-emitting radio neu- 
tron stars that the full HTRU survey will uncover. 
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Figure 1. A diagrammatic representation of tlie HTRU survey single-pulse pipeline. Each of up to 13 "single beam" processes associated 
with one multibeam pointing are run simultaneously, distributed onto several compute nodes and producing single-pulse and periodicity 
search results for that beam. The multibeam monitor tracks the status of the individual beams, collecting single-pulse search results, 
performing interference excision based on multiple-beam occurrence, and producing single-event diagnostic plots for each event in the 
pointing not tagged as interference. 



2 THE HTRU SURVEY'S SINGLE-PULSE (SP) 
SEARCH 

2.1 Data Analysis Pipeline 

Thus far all processing and analysis has been carried out us- 
ing the "Green Machine" supercomputer at the Swinburne 
Centre for Astrophysics and Supercomputing. The pipeline 
employs several standard search techniques (e. g. dedis- 
persion and box car matched filtering, mult i -beam coinci- 
dence matching; Cordes fc McLa ughlin 2 0031 : iDeneva et al.l 
|2009I :[k eane et al.ll2010l : BSB), and several novel techniques 
(e.g. independent-event identification and inspection). The 
basic methodology of the SP search is based on that de- 
scribed by BSB. 

There are two main bottlenecks in the HTRU survey's 
single-pulse search. The first is the computational load de- 
manded by the dedispersion of our data, and the second is 
the number of candidates which need to be manually as- 
sessed. The single-pulse data reduction pipeline design in- 
cludes schema that aim to abate both of these. First, we 
minimise the net computational time of the periodicity and 
SP searches by integrating the SP search into the hitrun 
processing pipeline introduced in HTRU Paper 1. By do- 
ing this, the interference excision and dedispersion described 
in HTRU Paper 1 is done only once, producing each data 
stream's results simultaneously for use in the single-pulse 
search and the periodicity search. As described below, the 
time series search for single pulses is performed while each 



time series is still in the computer's "random access mem- 
ory", and for each 9-minute data stream adds <2.5 minutes 
of processing time onto the total ~3-hour timescale for the 
HITRUN Fourier-based pipeline. We stress that the added 
computational load of including a single pulse search to a 
traditional Fourier search is negligible. 

We address the issue of large candidate numbers by 
identifying independent events across multiple dispersion, 
boxcar trials, and beams, and by performing a simple 
threshold-based interference zapping. The simultaneous pro- 
cessing and viewing of data from multiple beams decreases 
the occurrence of candidates caused by low-level, local inter- 
ference and greatly increases the speed at which the manual 
assessment of a candidate can be made. However, due to lim- 
ited computational resources, we cannot store all 13 beams 
on one node, requiring the use of a monitor to track and 
collect beams from each pointing as they are processed. 

The single-beam SP search and the multibeam data 
monitor are described below. Figure [T] gives a schematic dia- 
gram for the pipeline, showing one beam and the multibeam 
monitor. 



2.1.1 Single-beam Processing 

For each observation, a multibeam monitor is initiated and 
the available beams for that pointing are submitted in se- 
quence for processing on the supercomputer, each processed 
independently using the hitrun pipeline. For this pipeline 
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Table 1. HTRU survey observing parameters and an indication 
of the processed portion of the survey. Full survey parameters are 
given in HTRU Paper 1. Notes: ^The typical single-beam field of 
view is calculated within the average half-maximum beam width 
(~ 0.044 deg2). 





Latitude region 


Parameter 


Inter. 


High 


Avg. pointing duration (s) 


540 


270 


Total time processed 


3490 h 


140 h 


Total area processed^ 


7671 de{ 


5^ 591deg2 


Fraction of total survey 


23.5% 


0.39% 


Center frequency (MHz) 




1352 


Total bandwidth (B, MHz) 




340 


Typical usable Afghan 




870 


Sampling time (tsamp, 




64 



we begin with 2-bit filterbank data with the sampling prop- 
erties of which are summarised in Table[T] Pre-analysis radio 
frequency interference excision is done first, as detailed in 
HTRU Paper 1. In brief, this excision creates and uses both 
a frequency and a time mask that flags periodic-interference- 
affected channels and time samples with a DM = Opccm""^ 
signal above a signal-to-noise ratio (SNR) of more than 5. 
Flagged time samples are replaced in the filterbank file with 
noise drawn from adjacent unflagged samples, while offend- 
ing frequency channels are blanked. While the time-domain 
interference excision weakens the SNR of bright signals of 
DM < 0.12(^1/64 /is) pccm""^, where w is the pulse's width, 
it allows us to produce and inspect candidates down to DM 
~lpccm~^ while maintaining a manageable false detection 
rate at the manual inspection stage. 

We perform dedispersion at 1196 trial DMs over the 
range to 1000 pccm~^ using the dedisperse_ALL pro- 
gramQ After reading the input filterbank data into com- 
puter memory, time series are formed for each DM trial; the 
efficient dedispersion aspect of the dedisperse_all code is 
described in detail by HTRU Paper 1. As the time series 
data stream at each DM trial is produced, it is searched for 
single pulses with parameters as in BSB. Our input search 
parameters differ only in that the boxcar filter used in our 
search ranges in size A*';, = 1 to 512 samples, and that events 
separated by more than 10 samples (640 /is) were considered 
independent. All significant (as defined by our pre-set SNR 
threshold, mt > 6, detailed below) events in each boxcar 
and DM trial are recorded and stored in a database. At the 
completion of dedispersion of all trial DMs, dedisperse_all 
performs a temporal coincidence matchin g (similar to the 
"frien ds-of- friends" method devel o ped b y iHuchra fc Gelleil 
l|l982l ) and used bv lDeneva et all (|2009l l and BSB to iden- 
tify time-coincident event clusters). Each event is defined 
by the parameters (DM, boxcar filter, width, time) at which 
the signal-to- noise is found to be the greatest, and consists 



of members at other DM and boxcar trials found to be co- 
incident. 

The presence of man-made interference (i. e. radar com- 
munications, satellite and aircraft transmissions, on-site 
hardware, and a number of other sources) during observa- 
tions causes the noise in our time series data to follow a 
non-Gaussian distribution. These signals are typically either 
not dispersed, or dispersed to a level that is undetectable 
in our data, and therefore the zero-DM timeseries interfer- 
ence mitigation performed before processing removes the 
bulk of these signals. After interference pre-removal, how- 
ever, some low-level inteference signals remain, and we em- 
ploy pre-set search thresholds to balance the false detec- 
tion rate. For most data, the interference pre-removal al- 
lows time series at trials above DM = 1.5pccm~^, in the 
absence of our target astrophysical signals, to be roughly 
Gaussian-distributed. At lower dispersion trials, the false 
detection rate becomes unmanageable, and we therefore 
set a dispersion-based threshold to reject all events found 
with a peak SNR at DM < 1.5pccm~^. For each 9-minute 
intermediate latitude beam, we produce 1196 dedispersed 
time series of A^s — 8450000 samples. Our total number 
of searched data points for each file consists therefore of 
1196 X ELo(-^''/2') ^ 2.02 X 10^^ In Gaussian-distributed 
data, at our SNR threshold of mt = 6 we would expect 
approximately 50 independent random noise detections per 
file. We attempt to filter these events by automatically re- 
jecting events which have less than three associated mem- 
bers. This is an effective method of Gaussian event rejec- 
tion; noise peaks at a SNR of 12 or below are not likely 
to exhibit SNR > 6 in more than two trials, assuming the 
peak's signal between DM steps and boxcar trials drops by 
at least -\/2- Furthermore, the number of statistically random 
events above SNR = 12 is negligible, even when consider- 
ing analysis of the full HTRU intermediate latitude survey. 
This strategy impacts weakly on our sensitivity to events 
of w — Nbtsauip that have SNR < 12. However, in concert 
this event-match-based candidate rejection typically reduces 
the number of candidates by a further factor of ~100-1000, 
and much more for data badly affected by interference. In 
Fig. (2] we show the two visual inspection plots used for man- 
ual candidate descrimination. For the observation displayed, 
the telescope was pointed at PSR J1129-53, a known RRAT 
(BSB). The six pulses emitted by the pulsar in this point- 
ing were correctly discriminated by the software from zero- 
dispersion interference and from spurious peaks in the data. 

All non-interference events are then written to disk by 
DEDISPERSE_ALL, producing one ASCII file per event, in 
which all members of that event are listed. At this stage, 
the single-pulse events for the beam are gathered by the 
multibeam data monitor associated with the beam's point- 
ing. The time series are written to disk, stored for use in 
the periodicity search and for later access by the multibeam 
data monitor. 



^ This software and the C-l — h GTOOLS library, which is drawn 
from by dedisperse_ALL for single-pulse functionality and con- 
tains the SP search algorithms, candidate matching func- 
tions, and SP candidate type classes, are available from 
|http:/ /www.git hub.com/swinlegion| 



2.1.2 Multibeam Data Monitor 

As single-pulse events are produced for all the beams asso- 
ciated with one pointing, the multibeam monitor for that 
pointing collects the events and performs a temporal co- 
incidence match of identical form to that done by dedis- 
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Figure 2. Example diagrams used in the manual inspection stage, (a) FoUowing lCordes &: McLaughliiil | |2003|) , a single-beam summary 
diagram showing a position near PSR J1129— 53. 13 such plots are made for each pointing. Here events are coloured according to their 
auto-identified event type, as determined by the process described in t|2.1.1l Events flagged as zero-DM interference are black, those 
identified as gaussian peaks are red crosses, and candidate events are highlighted in blue/gray. The upper middle panel shows the SNR 
of each event plotted against DM, and the left panel shows a histogram of the net number of candidate (blue/gray) and all (black) 
detected event members in the pointing. The bottom panel shows the DM vs. time, with each member's SNR coded in the size of the 
plot point. In this example, six candidate pulses from PSR J1129— 53 have been correctly discriminated from the noise and interference 
in the observation. (b,c) Here we show single-event plots that correspond to the second and third pulses detected from PSR J1129-53 
in panel a. Panels show (clockwise from top left): the dedispersed timeseries in 13 beams; the SNR vs. DM (left subpanel) and boxcar 
filter trial (right subpanel) of the event; a false-colour image representing the signal power over frequency (MHz) and time since the 
beginning of the file (s); dedispersed time series at DM = Opccm"^ (bottom subpanel), and three DM trial steps around the brightest 
detec ted DM (upwards from bottom subpanel). The structural deviations in the SNR vs. DM curves in (b) and (c) from the predicted 
ICordes fc McLaughlini (^2003i') curves are the combined result of noise, and both time- and frequency-dependent pulse structure. 
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PERSE_ALL, however it is performed for events from different 
beams. 

The Parkes 20-cm multibeam receiving system allows 
the simultaneous observation of 13 positions on the sky with 
approximately 30 arcminutes between beam positions, and 
each beam has a sensitivity fall -off on scales of <30 arcmin - 
utes from its pointing centre l|Stavelev-Smith et al.lll99^ . 
Pointlike radio sources boresight to the telescope pointing 
direction will therefore typically be detected in a maximum 
of three beams of the receiver. Particularly luminous objects 
like the Vela Pulsar, which emits single pulses with peak 
flux densities of up to ~60 Jy, may be detected in up to 
~7 beams when ideally positioned in the multibeam field. 
Typically, signals of sufficient brightness to appear in all 
beams at similar intensity, i. e. through a far sidelobe of the 
telescope, are generated by terrestrial or near-Earth sources 
(such as satellites or aircraft). Therefore, as in BSB, we do 
not inspect candidates which were detected in more than 
nine of the thirteen beams. This filter decreases our sensi- 
tivity to the terrestrial "Perytons" of iBurke-Spolaor et al.l 
l|201lll . however there is significant benefit gained through 
decrease in the false detection rate per pointing. A search of 
the data aimed specifically at detecting Perytons is planned 
for the near future. All events occurring in nine beams or less 
are imaged as described in BSB, collecting the relevant in- 
terbank and time series data from the original on-disk loca- 
tion of the single beams. When both single-pulse and pulsar 
searches have been completed, the data is freed for removal 
from disk. 

Finally, a pointing's result plots (i. e. those for individ- 
ual beams and events as in Fig.[2)| are manually assessed to 
determine whether they contain a detection of interest. The 
"beam summary" plot has superiour sensitivity to objects 
emitting multiple pulses with signals at or just exceeding 
the detection threshold, while the single-event plots allow a 
user swift discrimination between interesting candidates and 
falsely-identified interference or noise. During the manual in- 
spection stage, results are also compared with the most up- 
to-date version of the Australia Telescope National Facility 
Pulsar Catalogue originally published bv iManchester et al.l 
l|2005[l . 

2.2 Search Sensitivity 

The sky coverage of the IfTRU survey includes regions cov- 
ered by previous surveys that have been searched for single- 
pulses. Particularly for the IfTRU intermediate latitude sur- 
vey, we have direct overlap with the areas searched by BSB 
and PKSMB. Ifere we calculate our sensitivity to transient 
events and make a comparison to these surveys. 

Instrumental and interstellar pulse broadening serves to 
weaken the SNR of the pulse from its signal at the pulse's 
intrinsic width, Wi. All instrumental broadening effects are 
in principle avoidable, therefore the back-end design of any 
observing system aims to reach the scatter- broadening limit, 
as this is the theoretical best that a pulsar/transients sur- 
vey can achieve. As noted in HTRU Paper 1, for a large 
range in DMs the improved time and frequency resolutions 
of the HTRU survey causes pulse broadening from inter- 
stellar propagation to dominate over intrachannel dispersion 
broadening (tch)- Usi ng the emp i rical D M-scatter broaden- 
ing relation given bv lBhat et al.l (|2004l ). scatter broadening 
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Figure 3. The HTRU survey's single-pulse search sensitivity to 
pulses of various intrinsic width, compared to the performance of 
the Parkes 20-cm analogue filterbank used by BSB (note: PKSMB 
had tsamp = 250 /iS whereas BSB had tsamp = 125 /is; thus 
the PKSMB curve Ues sUghly above BSB's). BPSR's higher fre- 
quency resolution most markedly improves our sensiti vity to nar- 
row p ulses at low DMs; above DM ~450pccm~^, the lBhat et al.l 
1120041 ) model for interstellar scattering dominates the instrumen- 
tal broadening for both surveys. Because of our limited search 
range in boxcar matched filter sizes, there is a sharp decrease in 
our sensitivity for pulses of duration > 32 ms. 



(tscat) typically dominates over hardware-induced smearing 
above a DM of ~230pccm~^ for the HTRU survey (note, 
however, that the scatter of individual pulsars in the Bhat 
et al. relation means that individual objects may deviate by 
up to an order of magnitude) . To calculate our sensitivity to 
single pulses, we take these pulse-scattering and broadening 

effects into account^ 

Following e..g:. iLorimer fc Kramej (|2005l ). the observed 
pulse width depends on Wi and various pulse-broadening ef- 
fects by w = {wf + + i?cat + *Lmp)°'^- Note that our step 
size between DM trials is chosen such that the broadening 
due to an error in DM is small compared to other broadening 
effects. If the observed peak flux density for a pulse is given 
by 5pcak, our sensitivity limit for the intrinsic peak flux den- 
sity of single pulses with observed duration w = A'^b x tsamp 
is then 



5'ilim ^ '5'pcak 

Wi 



'peak - , (1) 

GyJNpwB 

fc.f. iLorimer &: Kramerir2005l ) . where T^ys — 23 K is the sys- 
tem temperature of the multibeam system, G is the telescope 
gain (ranging from 0.735-0.581 K Jy~^, from the central to 
outer beams), and A^p = 2 is the number of summed po- 
larisations in the data. /? is a factor of order ~1 that is 
included to represent signal losses due to system imperfec- 
tions (for_^PSR^s_two2bit_digife /? ~ 1.07, following 
from lKouwenhoven fc Voutel200ll l. In Fig.El we indicate our 
sensitivity to single pulses of various durations as a function 
of DM, and in comparison to the sensitivities of PKSMB 
and BSB. 

The 2-bit digitisation levels of the BPSR instrument 
and large number of frequency channels result in the in- 
creased dynamic range capabilities over the previous Parkes 
analogue filterbank. We made empirical measurements of the 
mean (/i) and standard deviation (a) of HTRU survey data 
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judged by eye to not be strongly affected by interference. 
Using these values, we calculate the highest achievable SNR 
in the HTRU survey data for pulses with w = fsamp to be 
41.5. This is in good agreement with the SNR of saturated 
interference pulses observed in DM = Opccm"'^ time series. 
For the 1-bit, 96-channel analogue filterbank, the maximum 
SNR may be calculated analytically, using a binomial prob- 
ability distribution. The theoretical maximum SNR of this 
system at w — fsamp is then 9.6; again, this agrees with 
saturated interference signals observed in data from this 
system. The differences in dynamic range for these observ- 
ing backends do not strictly impact the flux sensitivity of 
a single-pulse search, however the HTRU survey's increase 
in dynamic range will affect the manual inspection stage by 
allowing detections of potentially very bright, narrow pulses 
to be more clearly discernible from noise or non-Gaussian 
statistics in the observation. This effect is more acute for 
pulses of width close to the sampling time of each survey. 

The improved time/frequency resolutions and dynamic 
range of the HTRU survey is clearly beneflcial for the detec- 
tion of the narrow single pulses of pulsars and RRATs, and 
extragalactic impulse emitters with durations w < 10 ms 
and DMs of DM < 360pccm"^. For DM > 360pccm"^, we 
do not expect to see a large number of new RRAT discoveries 
for several reasons; first, we do not have a considerable in- 
crease in sensitivity over the previous surveys for these DMs, 
and the number of known RRATs above this DM repre- 
sents only ~10% of the known population. Additionally, for 
l&l < 5, the 35-minute pointing duration of PKSMB searches 
afford a greater probability of detection of low-pulsation rate 
events over our 8.5-minute observations of these areas, and 
we offer minimal improvements for the capability of detec- 
tion of these objects. At higher latitudes, for the bulk of 
pointing positions in the HTRU survey the number density 
of electrons in the Gala:xy is in sufficient to produce a DM 
much greater than 360 pc cm" l|Cordes fc Lazidl2002l ). For 
these reasons, we expect most new rotating radio transient 
and high nulling-fraction pulsar discoveries to be at DMs 
DM < 360pccm"^. As detailed in 33] and gT] this is true 
for the initial discoveries. 



2.3 Current Processing Status 

The current status of the HTRU survey's processing is sum- 
marised in the upper panel of Table [T] The processing re- 
ported in this paper includes only pointings from the HTRU 
intermediate and high latitude surveys. The survey coverage 
to date is not contiguous over all Galactic regions. We in- 
dicate the total observation time of the processed pointings 
as a function of sky position in Figure [4] In total, 22,273 
and 1,718 beams have been processed out of the 95,056 and 
443,287 total beams of the intermediate and high latitude 
surveys, respectively. Approximately 160 incomplete survey 
pointings of shorter duration have also been processed. As 
implied by Fig.|4l the bulk of observing time reported in this 
paper has been spent at Galactic latitudes |6| < 15. 

2.4 Candidate Ranking and Follow-up Strategy 

We employ a follow-up prioritisation scheme for our single- 
pulse candidates to provide maximum scientific return with 



available observing time. Based on the strength of the candi- 
date in various panels of the manual inspection plots and the 
number of detections at a similar DM in a pointing (Fig.[21), 
candidates are ranked as multi-pulse (two or more strong 
pulses detected at similar DM), strong (^2 weak pulses, 
or one convincing pulse detected, that is with a SNR > 7 
and exhibiting a characteri stic signature in the SNR vs. DM 
plot, e.g. as explored bv ^ Cordes &: McLaughlinI |2003| ). or 
weak (a DM > Opccm"'^ detection appearing pointlike 
in the multibeam field, but poorly supported by other in- 
spection panels and/or has SNR < 7). Weak candidates 
outnumber those ranked strong and multi-pulse by a fac- 
tor of more than two. All candidates that are not found 
in available archival data with pointings near the discovery 
position are reobserved for a duration equal to that of the 
discovery pointing. If weak detections are not seen in the 
initial follow-up pointing, they are not tracked in further 
follow-up observations. We do not include these detections 
in this report due to their highly uncertain nature, however 
we preserve them in a database for posterity, considering 
that borderline detections in future surveys of the same sky 
area could cross-correlate their findings with our weak, un- 
confirmed detections. This will be possible with the HTRU 
low-latitude survey, and potentially with future Australian 
Square Kilometre Array Pathfinder/Square Kilometre Ar- 
ray surveys. Multi-pulse and strong candidates not seen in 
initial follow-up are observed again for a duration of three 
times the discovery pointing. Our ranking levels are set such 
that the strong and multi-pulse categories are highly unlikely 
to contain detections of spurious noise. We include all such 
candidates here, with objects not detected in follow-up ob- 
servations distinguished in Table [2] by an upper limit on the 
objects pulsation rate (x, quoted as number of pulses per 
hour); this reflects that non-detections of these discoveries 
are used only to place limits on x for the object. While the 
limits are well within the range of pulsation rates for known 
RRATs, these objects essentially remain unconfirmed, and 
care must be taken in including these objects in statistical 
or categorical transients studies (note, however, that PSR 
J 1424-56 is conflrmed by the reported detection of PSR 
J1423-56 by Keane (2010b), as discussed below. 



3 EARLY DISCOVERIES &z DETECTIONS 
3.1 New Discoveries 

Table [2] lists basic parameters of the 13 new and conflrmed 
objects discovered thus far using the HTRU survey's single- 
pulse pipeline. There are several objects worthy of speciflc 
note. 

PSR J0410-31, at Galactic coordinates I = 230.6°, 6 = 
—46.7°, is the only conflrmed single-pulse discovery arising 
thus far from the high-latitude survey. The DM of this pulsar 
is low, with DM = 9.2pccm~^, lower than 98.9% of known 
radio pulsars. Estimati ng distance to the new discoveries us- 
ing their DMs and the ICordes fc Laiol (2002) electron den- 
sity model for the Galaxy (hereafter NE2001), we flnd that 
this is also the most nearby object discovered in our search 
(d = 0.51 kpc). 

PSR J0912— 38 was also detected in the HTRU sur- 
vey's Fourier-domain search. Bright single pulses were only 
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Figure 4. A Galactic map showing new HTRU survey single-pulse discoveries, and observing time spent per (, 6 cell. The stars indicate 
detections from which multiple pulses have been detected, and triangles indicate sources which have only yet exhibited one strong pulse. 
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Table 2. Properties of the objects discovered in this search. Columns: (1) Name based on J2000 coordinate indicates a candidate from 
which only one clear pulse detection has yet been made. Such objects must be interpreted with care; see notes on these and the objects 
with limits only on x in HMJ. * PSRs J1307-67 and J1423-56 are believed to be confirmations of PSRs J1308-67 and PSR J1424-56, 
respectively, of Keane (2010b); see notes about these object in Section|3l (2,3) J2000 right ascension and declination of pointing centre 
for detected beam. Only for J1854— 1557, the position is as derived from the pulsar's timing solution; (4) the best-fit period, where 
measurable, with the error on the last digit in parentheses; (5) pulsation rate x = ^ph~^; (6) best-fit DM and error in pccm"'^; (7) 
observed pulse width at half-maximum of the brightest pulse; (8) peak flux density of the brightest detected pulse (calculated using 
Speak in Eq.[lJ. 



PSRJ 


Pointi 
RAJ 


ng centre 
DECJ 


P (s) 


X 

(h-i) 


DM 


(ms) 


^pcak 

(mJy) 


J0410-31 


04:10:39 


-31:07:29 


1.8785(2) 


107 


9.2(3) 


18 


470 


J0837-24 


08:37:44 


-24:47:48 




5 


142.8(5) 


1 


420 


J0912-38 


09:12:27 


-38:48:34 


1.5262(1) 


32 


73.3(5) 


6 


190 


J1014-48 


10:14:18 


-48:49:42 


1.5088(2) 


16 


87(7) 


21 


140 


J1135-49t 


11:35:56 


-49:25:31 




^1.3 


114(20) 


9 


120 


J 1216-50 


12:16:20 


-50:27:01 


6.355(9) 


13 


110(20) 


9 


130 


J1307-67* 


13:07:41 


-67:03:27 


3.65120(8) 


11 


47(15) 


32 


70 


J1424-56* 


14:24:23 


-56:40:47 




s:7 


27(5) 


22 


110 


J1541-42t 


15:41:55 


-42:18:50 




s;7 


60(10) 


4 


150 


J1549-57 


15:49:05 


-57:21:37 


0.7375(3) 


73 


17.7(3.5) 


4 


210 


J1709-43t 


17:09:47 


-43:54:43 




s;7 


228(20) 


3 


240 


J1854-1557 


18:54:53.6 


-15:57:22(20) 


3.4532(1) 


25 


160(25) 


65 


50 


J1925-16 


19:25:06 


-16:01:00 


3.8858(2) 


^6.5 


88(20) 


10 


160 



detectable above the noise for <30 seconds in both the dis- 
covery and follow-up observations. When the data are folded 
over pulses not detected by the single-event search, the pul- 
sar remains detectable with a signal-to-noise ratio of >5; this 
behaviour is not consistent with scintillation. The abrupt 
change in flux density and the occurrence of these short- 
duration flares in both observations lead us to suggest they 
are intrinsic to the star. 

The pulses of PSR J1014-48 were only detected in one 
cluster spread over 16 rotations of the neutron star. Like 
PSR J 1825-33 of BSB, no outbursts have yet been detected 
in follow-up observations. 

We are confldent that PSR J1307-67 is the PSR J1308- 
67 of Keane (2010b), who reports a solitary pulse discovery 
at DM = 44± 2pccm~'' and a position with errors in agree- 
ment with our pointing position. Upon searching the archival 
data closest in position to our discovery, we found 7 pulses 
with a signal-to-noise ratio of >5, from which the period re- 
ported in Table[5]was determined. Similarly, we are confident 
that PSR J1424-56 is the PSR J1423-56 of Keane (2010b), 
who reports DM = 32.9 ± l.lpc cm"^ and x = 3.4 h"\ both 
consistent with our findings. Timing solutions will solve the 
positional and naming ambiguities of these objects. 

Broad-bandwidth amplitude modulation features are 
seen in the frequency-dependent structure of the single 
pulses of PSR J 1549-57 (see Fig.[5l), the bandwidth of which 
are inconsistent with the predicted scintillation bandwidth 
of NE2001 by four orders of magnitude. The origin of this 
feature is unknown, however despite the significant disagree- 
ment with the NE2001 prediction, it is likely to be scintilla- 
tion if the detected emission is produced by a neutron star. 

The predicted scatter broadening of PSR J1709-43 
(130 ms at 1.375 GHz for the NE2001 model) is much greater 
than the half-maximum pulse width detected here. The sin- 
gle pulse from this object had a signal to noise ratio of 
10, a visible dispersion curve , a signal-to-noise versus DM 
curve well-fit to the model of lCordes fc McLaughTinI (|2003l ) 



IX 




0.55 



Figure 5. PSR J1549— 57: (a) A pulse stack of 16 single pulses, 
(b) A dedispersed waterfall image (lower panel) with dedispersed 
and band-integrated pulse profile (upper panel); these result after 
folding only over the 16 stellar rotations shown in (a). 



for a genuinely dispersed pulse, and exists in a relatively 
interference-free observation. These conflicting properties 
indicate that either the impulse's origin is not celestial, or 
that the NE2001 scattering model for this Galactic position 
is incorrect; the latter is highly possible to be the case here, 
as scattering is highly dependent on line-of-sight. A redetec- 
tion of this candidate would conclusively clarify the nature 
of this candidate, however follow-up has not yet been per- 
formed. 

PSR J1854-1557 appears to be exhibiting periodic 
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Figure 6. A pulse stack showing single pulses of the periodically nulling pulsar J1854-1557. 



nulling, drifting, and mode-changing behaviours (Fig.|6]). 
Despite a high nulling fraction, this object's relatively fre- 
quent emission renders it detectable in a Fourier search with 
a signal to noise ratio of ~17 in the 8.5-minute survey and 
follow-up pointings. This object was also detected in the 
Fourier pipeline, as reported in HTRU Paper 1. 



3.2 Redetections of Known Pulsars and RRATs 

The single-pulse pipeline has detected approximately 55% of 
the known pulsars that could be detected through Fourier 
searches of the HTRU survey data. This is an unexpectedly 
large fraction, and it offers improvements over the single- 
pulse detection rates presented in BSB and PKSMB, per- 
haps a testament to the increased sensitivity to faint and 
narrow single pulses that this survey's hardware affords. The 
single-pulse properties of these objects will be analysed in 
future work, so we will not detail them here. 

We have processed 24 observations with point- 
ing positions within 6 arcminutes (our approximate 
half-powe r beam width) of published RRAT positions 
fi. e. from iMcLaughlin et al ll2006l. iMcLaughhn et al.|[2009l . 



iKeane et al.ll20ld . iBurke-Spolaor fc Bailedi2010l ). Of these 
24, 13 RRATs were redetected. For 22 of 24 observations, 
the number of pulses seen from each object was consistent 
with the expected number based on previously published 
pulsation rates in the discovery observations of these objects 
(again, as drawn from the BSB and PKSMB publications) 
and the observation length in our data. The two notably 
different detection rates were of RRATs J1753-38 (BSB; 4 
expec ted and > 15 seen) and J1819-1458 (jMcLaughlin et al] 
l2006l : 3 expected and 10 seen). PSR J1819-1458 is known 
to exhibit small variatio ns in pulsation ra te (and larger ones 
associated with glitches. iLvne et al.ll2009l ). and the detection 
rat e observed here is consistent with the range in x reported 
bv lLvne et all l|2009l ). For PSR J1753-38, the higher detec- 
tion rate, detectability in a Fourier search with a signal-to- 
noise ratio of ~13 (whereas it was undetectable in a Fourier 
search of its original discovery data), and the slight increase 
in Smax exhibited in our observations indicate that: 1) this 
pulsar's emission appears to simply be highly modulated 
rather than nulling, and 2) our observing position of the 
source is improved from that reported by BSB. Our posi- 
tional offset was 0.1 degrees, suggesting an improved posi- 



tion of right ascension and declination 17:53:09, -38:52:13 
(J2000). 



4 DISCUSSION 

4.1 Our Discoveries in Archival Pulsar Surveys 

All of our discoveries except for PSR J0410-31 and PSR 
J0837-24 lie in regio n s pre vi ously covered b y the sur- 
veys of Jacobv et al.l l|2009l) : lEdwards" et all I2OOII ) and 
[Manchester et al. I l|200lh . The presence of pulses in archival 
data can add valuable information about pulse rate changes 
over time, and potentially add data points to their timing if 
a phase-coherent timing solution is obtained. The detection 
of J1307-67 and J1424-56 in archival data have already been 
discussed in Sec. 13.11 We inspected data in the archival sur- 
veys within one half-power beam width of our other objects' 
positions to determine the detectability of our discoveries in 
the archival data. None of our objects were detected in these 
pointings in single pulses. However, while PSR J 1854-1557 
is not detectable in a Fourier search, it is marginally de- 
tectable when the archival data are folded over the period 
given in Table [21 

We investigate the reason for the non-detection of the 
remaining objects by first estimating the signal loss due to 
the use of the analogue multibeam filterbank system used 
by BSB and PKSMB, calculating the SNR that a pulse of 
equivalent width, DM, and 5pcak (as listed in Table[2| would 
have shown in the archival data. This analysis suggests that 
even the brightest detected pulses in our data would result in 
only a marginal (SNR < 6) detection in the analogue filter- 
bank data for PSRs J1135-49, J1541-42, and J1709-43. For 
two out of four of the remaining objects at |6| > 5°, for which 
the archival data length was 4.4 minutes, the non-detections 
are accounted for by the low probability that a pulse would 
occur during the time span of the observation. The remain- 
ing |6| > 5° objects, PSR J0912-38 and PSR J1014-48, had 
clusters of sequential pulses rather than a smooth distribu- 
tion of single pulses. Although the per-hour pulsation rate 
is high for these objects, the duration of on-activity is short, 
and the spacing between sequential pulse outbursts is >8 
minutes. This indicates that the non-detection of pulses is 
consistent with a decreased probability of the occurrence of 
a pulse outburst during the archival observations. The ef- 
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feet of on- and off-timescales on object detectability will be 
discussed further in Sec. 14.21 

The only remaining undetected object is PSR J1549-57. 
We have not yet performed follow-up observations on this 
object, and so the archival data, in which we would have 
expected to see detections of SNR > 7, serves only to place 
limits on the pulsation rate and duration of off-activity in 
these objects. The non-detection in two archival pointings 
places a strong limit (zero pulses in ^70 m) on these values. 



4.2 Activity Timescales in Sparsely-emitting 
Neutron Stars 

Approximately 1/3 of our single-pulse pipeline discoveries 
exhibit distinct periods of on-activity, marked by bright se- 
quential pulses that are separated by longer intervals of ei- 
ther genuine nulls, or a decrease in intensity level sufficiently 
large that the object is undetectable through single pulse 
or periodicity search techniques during the null state. The 
pulsation rate, Xi that is typically quoted for single-pulse 
discoveries poorly represents the emissivity of such objects. 
It provides insufficient information to assess either the prob- 
ability of discovering such a neutron star in a given observa- 
tion, or of detecting such an object during a reobservation. 

Considering this and the significant capability of the 
IfTRU survey to discover new genuinely sparsely-emitting 
objects (see i]4.4|) . we seek to more accurately incorporate 
our discoveries and RRAT populations in general into the 
range of nulling and emissivity timescales exhibited by ra- 
dio pulsars. We explore a relevant "intermittence parameter 
space" here, and review in H4.3I what effect survey param- 
eters have on the selection of various populations in this 
space. 

We note that the on- and off-stat es of nulling pulsars 
typically show characteristic l e ngths ( Herfindal fc Rankin! 
l2007l : iRedman fc RankinI |2009| : IWang et al.ll2007l ). paranT 
eterised below by ton and toS, respectively. This description 
of neutron star intermittence is more physically representa- 
tive for pulsars of various nulling fractions, better reflects 
the timescales associated with possible windowing phenom- 
ena, and allows a more accurate exploration of selection ef- 
fects for surveys of various length (see H4.3p . Here we use 
a definition of tpfi that is similar to the "null length" of 
IWang et al.l ((2003). We define toff as the average time be- 
tween the first pulse whose signal drops below a set thresh- 
old (we use SNR > 5) and the next above that threshold. 
We define ton as the average number of pulses (Non) above 
the same threshold, times the rotational period of the ob- 
ject. Obviously, such values are only valid for genuine nulling 
pulsars, that is this analysis is inappropriate for highly mod- 
ulate d pulsars with lo ng-tailed pulse energy distributions 
(e.g. IWeltevrede et aU 2006). Uninterrupted data spans of 
T S> ton + tos containing single-pulse detections of high sig- 
nificance (that is, where the single pulse energy distribution 
of the on-pulses is such that all pulses exceed the thresh- 
old) would provide the most accurate measurements of these 
quantities; for our discoveries, most observations provide 
robust single-pulse detections, although we are insensitive 
to timescales ^ 9 minutes. We nevertheless make estimates 
of (or place limits on, in observations where no timescales 
are measureable because they appear greater than our data 



span) the average ton, toft for our discoveries based on all 
available data. 

The nulling fraction of a pulsar is only valid in this 
analysis when measured over a timescale far exceeding one 
activity cycle of the object, which is given by a length 
Ta = ton + tod- Below we generically represent the nulling 
fraction measured over "infinite time" as foo, where ac- 
cordingly, the pulsar's on fraction measured over a time 
r ^ Ta is given by ^ = 1 — fao- As an example, the pul- 
sar PSR B1931+24 exhibits periods of activity lasting ~5- 7 
days, separated by nulls of ~30 days (|Kramer et al.|[2006l ). 
While the nulling fraction measured over short observational 
timescales is therefore typically either 0.0 or 1.0, we find that 
/oo 0.85. 

In Figure [71 we plot ^ as a function of toft for the 
single-pul se discoveries from HTRU, PKSMB, BSB, and 
lDeneva"et a l. (200^ We al so plot all nulling pulsars anal- 
vsed bv IWang et al] (|2007l ) that have a non-zero nulling 
fraction. The Wang et al. sample was selected from the 
periodicity-discovered pulsars of the same survey as the 
PKSMB searches, and has poor sensitivity to nulls of ^ 30 s 
due to the averaging performed over 10-30 s timescales. The 
dashed grey lines in this figure show contours of constant 
ton. For a pulsar of period P, it is not straight-forward to 
identify on- or off-activity timescales of duration t < P; 
the shaded region below ton = 100 ms indicates the area for 
which ton is shorter than the period range of single-pulse 
search discoveries that emit at A'^on ~ 1. 

One critical point illust rated by this plot i s the contin- 
uous transition between the I Wang et al. 1 (|2007l ) nulling pul- 
sars and the populatio n discovered by single-pulse searches. 
[McLaughlin fc Corded I^OOS ) derive that neutron stars of 
period P will be more readily discovered in a single-pulse 
search, (i.e. having SNR ratio r = msp/mppT > 1) when 
C < 2^/P/T, leading Keane (2010a) to note that this ef- 
fect causes ^- and P-based selection effects for single pulse 
searches. For the HTRU intermediate latitude survey point- 
ings, pulsars with periods P < 6s should therefore be dis- 
covered more readily in the Fourier pipeline at ^ > 0.2, 
and thus it is not surprising that no HTRU intermediate 
latitude survey points lie above this line in Fig. [7] For 
the r = 35 min pointin gs of the Parkes Multibeam survey 
([Manchester et al.|[200ll ). r > 1 when ^ < 0.02 to 0.11, con- 
sidering a period range 200 ms to 6 s. Appropriately, this ap- 
pears to be the region in Fig.[7|where the nulling pulsars give 
way to the single-pulse discoveries. It is clear here that the 
on and off timescales, and the nulling fractions of the radio 
pulsars discovered in single-pulse searches exhibit a smooth 
distribution over many orders of magnitude with no obvi- 
ous distinction of the RRAT population represented by the 
single-pulse discoveries. A rigorous nulling fraction analysis 
for bright pulsars in the HTRU intermediate latitude survey 
is planned for the near future, and is capable of addressing 
whether the objects discovered by single-pulse searches rep- 
resent a smooth extension in the probability density function 
of pulsar nulling fractions. 

4.3 Detectability of Intermittent Neutron Star 
Populations 

The distribution of objects in Figure [71 is influenced by a 
number of selection effects, which we explore here both to 



12 S. Burke-Spolaor et al. 



0.001 



0.0001 





1 1 1 i| 1 




L=50% : 


■ X 

A 


A 




O 
L=95% ■ 


■ ® 

: VV ® ® 

■ % \ 


© 


A 


N ^ X BSB 

A PKSMB 

Denevaetal. (2009) - 
■ Reported here 
* Nulling pulsars 
O B1931+24 


\ 

e 

I 

z> 

DC 
1- 

. . 


oa 


A^ 

o 
z> 

DC 
1- 
,^ 


A \iJ> \^ \(J> 
1 " 1 \ 1 



10 100 1000 10000 100000 1e+06 1e+07 



toff (seconds) 



Figur e 7. Here we show the on fraction, 5, against the mean null length for single-pulse search discoveries and the pulsars of lWane et alj 
(|2007|). All elements highlighted in maroon indicate relevance to the HTRU intermediate latitude survey; the squares show single-pulse 
pipeline discoveries, while maroon circles indicate previously-known objects which were redetected in the HTRU intermediate latitude 
survey data ( i|3.2l l. Points with arrows designate an object for which ton has been measured but only a lower limit on the null timescale is 
known, and points with flat-headed arrows represent objects for which only one pulse has been detected (we have assumed P = 8.5 s for 
these objects, therefore the plotted point is an upper limit to g and a lower limit to toft)- The solid green lines indicate lines of constant 
nulling fraction as measured over infinite time (/oo, see i|4.2l l. and the dashed gray lines show constant contours of ton- The vertical 
dotted lines show the single-pointing observation length for the survey of PKSMB, and for the HTRU high (hi), intermediate (med), and 
low (lo) latitude surveys (see i|4.3ll . The shaded region below ton = 100 ms indicates the area for which ton is shorter t han the period 
range of single-pulse search discoveries, and is therefore ill-defined. B1931-I-24 is the "intermittent pulsar" published by iKramer et aU 
(|200y), and discussed in Sec. 14.31 



investigate the underlying distribution of activity cycles in 
neutron stars, and to determine the population to which the 
three HTRU survey components will be sensitive. 

Considering the general detectability of pulsars in this 
phcise space, the probability that at least one pulse will be 
emitted by a bright pulsar during an observation of length 
T is given by: 

r for = 1 

-fern — < 1 for T > toff (2) 

i (T - 5 T + etoff ) /i off for T < toff 

Note, however, that this probability is generally representa- 
tive and will hold for a sample of many sources, but will not 
always be strictly true for an individual source as the distri- 
bution of ton and toff for a single object is not a delta func- 
tion. Following this equation, the vertical lines in Fig.[7]mark 
the null length below which a sufficiently bright emitter has 
a probability of unity that it will be detected in the Parkes 



Multibeam survey, the HTRU intermediate latitude survey, 
and the (recently commenced) low-latitude survey, respec- 
tively. At null lengths greater than the observing time for 
these surveys, the probability of detection quickly decreases. 
Accordingly, given a population of pulsars with a flat under- 
lying distribution in toff or log(toff ), we would expect to have 
a higher discovery /detection rate of objects with toff < T. 
Focussing on the objects detected by the HTRU interme- 
diate latitude survey (highlighted in maroon in Fig.O, the 
distribution is roughly consistent with a flat distribution in 
log(tofr), with a ratio of 7:13 for objects with toff ^ 9min to 
objects with toff < 9min. We point out that the object of 
the highest toff discovered in each of the PKSMB, BSB, and 
HTRU intermediate latitude surveys is roughly a factor of 
three greater than the corresponding survey's length. This 
corresponds to Pcm ^ 0.33 at all nulling fractions. 

Equation [5] indicates that on average for toff < T, 
the probability a source would be detected is one, and for 
toff > T, it is directly proportional to ^. It is apparent in 
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Fig. [7| that there is a drop in Parkes Muhibeam detections 
in the region 100 s < toff < 35min, ton > 100 s. That is, 
if the underlying neutron star population has a flat distri- 
bution in both ton and tas or log(ton) and log(toff), there 
appears to be a deficit of objects in this region. For ob- 
jects with C ^ 0.1 (that is, those theoretically detectable in 
a Fourier-search), there is no explicitly stated reason that 
IWang et al. I ([2007) should have selected against this range 
of timescales. The paucity of objects occupying this region 
could be the combined result of a number of possible ef- 
fects. Nulling behaviour may have not yet been either ob- 
served or recognised, however we would expect this to only 
be the case for objects with ton 35 min, which encompasses 
only a small fraction of the region in question. It is possible 
that some low or intermediate nuUing-fraction objects with 
toff — 35 were outright rejected from the survey due to a 
non-appearance in follow-up observations. We consider this 
an unlikely cause, however, as the 6-minute grid pointings 
and subsequent 35-minute follow-up poi ntings (in case of 
a no n-detection in the grid observation, [Manchester et ahl 
I2OOII ) would have provided a sufficiently high redetection 
probability, particularly for objects with toff <^ T. We note 
that for objects with toff close to 35 minutes, the survey's 
follow-up observations might contain only a few pulses from 
the objects, which are detectable in single pulses, rather than 
in the Fourier search used to assess the follow-up pointings 
in the survey. However, this effect still does not account for 
the notable detection rate drop in the region in question. 

It is furthermore possible that the lack of objects in 
this region reflects an actual drop in the neutron star popu- 
lation for objects with ton ^ 300 s. This supposition is ques- 
tioned, however, by the exi stence of the " i nterm ittent pul- 
sar" B1931+24 reported bv lKramer et al.l (|2006t ). This ob- 
ject is isolated in this space, although it has been indicated 
that a handful of fur ther unpubl i shed objects oc cupy a sim- 
ilar region of Fig.E] (|Keithll2007l : lEatoughl [2OO9I ) . The exis- 
tence of these objects suggests that either there is a bridging 
distribution of neutron stars which we have not yet identi- 
fled due to an unrecognised effect of selection, or that in fact 
the intermittence of these objects is produced by a physi- 
cally distinct mechanism. Immense observing timescales are 
required for a survey to have a high probability of detect- 
ing a population between the ton ^ 300 s objects and PSR 
B1931-f24-like pulsars. Because data storage and analysis 
becomes problematic for pulsar surveys of longer timescale, 
we suggest that the brightest pulsars in this population may 
be more readily uncovered by the "transient-imaging" sur- 
veys that are being performed and planned on arrays; these 
surveys can sample timescales down to the correlator aver- 
aging time (typically >10s) and up to several hours to days. 
However, for pulsars with ton less than these timescales, pul- 
sar surveys remain the most robust probes of intermittent 
populations. 

4.4 Discovery Forecast for the Full HTRU Survey 

The rate of new single-pulse discoveries in the HTRU survey, 
particularly those in overlapping regions of previous surveys, 
hint at the potential for the full HTRU survey to uncover 
many new examples of sparsely-emitting neutron stars. For 
the intermediate latitude survey, we have detected a total 
of 26 new and known single-pulse emitters, suggesting de- 



tection rates on the order of 3 x 10~^ deg~^ at these lati- 
tudes. We have processed roughly 23.5% of the intermediate 
latitude survey to date, and if our discovery rate contin- 
ues, the full HTRU intermediate latitude survey data should 
contribute an additional estimated ~50 new discoveries of 
transient neutron stars. We note that the sky distribution of 
neutron stars should have a higher density at lower Galactic 
latitudes, and as seen in Figure[4]we have not yet processed 
most of the pointings at latitudes closest to 6 = 0. At higher 
latitudes we have only processed a small fraction of the sur- 
vey, and while tempting to take an extrapolated value as 
an upper limit, we note that our single discovery in only 
0.39% of the high latitude survey implies a sky density of 
~2 X 10~^deg~^, higher than expected when compared to 
the sky density at lower latitudes. Therefore, the detections 
are expected to be less than 250. 

Finally, the explorations of Sections 14.21 and 14.31 have 
accented the timescales of neutron star intermittence ac- 
cessable by the HTRU survey. In particular, the HTRU 
low-latitude survey (with its 70-minute pointing length) 
will mark a signiflcant increase in detection rate over the 
PKSMB survey for objects with 35 < toff < 70 min, and 
may be realistically expected to discover objects with null 
lengths up to toff = 3.5 h, if such objects exist (as deduced 
in i]4.3l in which it was noted that transients surveys have 
all been able to discover objects with toff values of up to 3 
times the survey pointing length; see also Fig. [7]). The low 
latitude survey will also provide data that may be used to 
investigate the possible deficit of low- nulling- fraction objects 
with ton > 300 s. Care must be taken during both the single- 
pulse and periodicity searches, analysis, and follow-up, such 
that if these objects do exist, they are not selected against. 
This will involve, for one, appropriate follow-up monitoring 
of the stronger Fourier candidates that are undetected in 
initial follow-up pointings. 



5 SUMMARY AND CONCLUSIONS 

We have presented the methods of and initial discoveries 
for the single-pulse analysis of the High Time Resolution 
Universe Survey. We outlined the design of the HTRU sur- 
vey's single-pulse pipeline, which functions efficiently along- 
side the Fourier-analysis pipeline; the single-pulse analysis 
furthermore employs a "friend of friends" single-event recog- 
nition algorithm and performs automated interference re- 
jection based on the dispersive and multi-beam signature 
of single events. The new digital backend used to collect 
HTRU survey data has afforded a factor of up to ~5 times 
improvement in sensitivitiy over previous surveys in over- 
lapping Galactic regions, and offers the most significant 
improvements for short (sub-ms duration) pulses at low 
(DM < 360pccm~^) DMs. 

Analysis of 23.5% and 0.39% of HTRU intermedi- 
ate and high Galactic latitude survey data, respectively, 
has resulted in the discovery of 12 and 1 new neu- 
tron stars. Much of the survey pointings covered Galac- 
tic regions tha t were previously survey e d and searched for 
single pulses jMcLaughlin et al.] 12009 : iKeane et all I2OI0I : 
iBurke-Spolaor fc Bailesll2O10D : 11 of our new discoveries lie 
within these regions, and for the nine of these that were not 
visible in the archival data, their non-detection was consis- 
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tent with either the signal degradation due to the use of the 
previous wider-channel analogue backend, or with the im- 
probability of a pulse being emitted in the archival survey 
due to a long null cycle. 

Finally, wc investigated the distribution of nulling and 
emissivity timescales for the new single-pulse neutron star 
discoveries and redetections in the HTRU survey data, and 
for the RRAT/nuUing population in general. We found 
that periodicity-discovered nulling pulsars and single-pulse 
search discoveries exhibit a continuous distribution across 
null/activity timescales and nulling fractions, building on ev- 
idence that many RRATs represent a tail of extreme-nulling 
pulsars. We found that there is an apparent decrease in ob- 
jects with emissivity cycles longer than ~300 seconds at in- 
termediate and low nulling fractions which is not readily ex- 
plained by selection effects, and note that the HTRU deep 
low-latitude survey will be capable of exploring whether this 
deficit is natural or an effect of selection. Lastly, we esti- 
mated that the full HTRU survey is capable of more than 
doubling the known extreme-nulling pulsar population, and 
will explore the neutron star population with nulling frac- 
tions exceeding 99.99%, and null lengths lasting up to 3-4 
hours. 
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